Orbital-Ordering Induces Structural Phase Transition and the Resistivity Anomaly in 

Iron Pnictides 
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We attribute the structural phase transition (SPT) in the parent compounds of the iron pnictides 
to orbital ordering. Due to the anisotropy of the dxz and dy^ orbitals in the xy plane, a ferro-orbital 
ordering makes the orthorhombic structure more energetically favorable, thus inducing the SPT. In 
this orbital-ordered system, the sites with orbitals that do not order have higher energies. Scattering 
of the itinerant electrons by these localized two-level systems causes a resistivity anomaly upon the 
onset of the SPT. The proposed orbital ordering also leads to the stripe-like anti-ferromagnetism 
and anisotropy of the magnetic exchanges. This model is quantitatively consistent with available 
experimental observations. 

PACS numbers: 61.50.Ah, 71.70.Ej, 72.15.Qm 



I. INTRODUCTION 

The structural phase transition (SPT) from tetragonal 
to orthorhombic symmetry around 150 is a ubiqui- 
tous feature in the parent compounds of the iron-based 
superconductors. Coincident with this transition is a re- 
sistivity anomaly (RA)^ in which the resistivity turns up 
slightly before a sharp drop at exactly the onset temper- 
ature of the SPT, TspT- For the 1 Ill-family, at a lower 
temperature, TsdWj a stripe- like anti- ferromagnetic spin 
density wave (SDW) formsi on the distorted lattice of Fe 
atoms, with the spins being parallel along the shorter axis 
and anti-parallel along the longer axis. However, for the 
122-family, the SDW develops at the same temperature 
as does the SPT, Tgow = Tspt^. In the 122-family^, a 
single first-order transition obtains instead of two sepa- 
rate second-order transitions in the 1111-family. Upon 
doping, superconductivity (SC) occurs leading to a ces- 
sation of the SPT, RA and SDW^i. Hence, all of these 
three phenomena should be closely related and share a 
universal mechanism. However, most theoretical work 
only focuses on the connections between the SDW and 
SC. The importance of the SPT and RA is somehow un- 
derestimated. The main objective of this paper is to 
explain the origin of the SPT and RA. 
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is that the SPT is driven by 



the onset of the stripe-like anti-ferromagnetism. Both 
first principles calculations^!^ and Landau-Ginzburg 
modelingsiSiiii have been used in this context. The fact 
that the two transitions are decoupled in the 1111-family 
is a limitation of this approach. Further, since the ori- 
gin of the SPT in their scenario is spin based, the onset 
temperature should be sensitive to an external magnetic 
field. However, experiments have shown that varying the 
magnetic field leads to no change in the onset tempera- 
ture of the SPT2. 

In this paper, we develop a microscopic theory of the 
SPT without involving the spin degrees of freedom. On 
our account, uneven occupations of the dxz and dyz or- 
bitals make the orthorhombic crystal structure more en- 
ergetically favorable, thus inducing the SPT. The opera- 



tive mechanism driving this ferro-orbital-ordering transi- 
tion is the lifting of the degeneracy between the dxz and 
dyz orbitals by the inter-site Coulomb repulsions. How- 
ever, it should be noted that other important factors, 
such as spin-orbit interactions^^ and couplings to the 
displacements of ligand atoms (As), also contribute to 
this process. In fact, spin-orbit physics appears to lie at 
the heart of orbital ordering in the manganites^. While 
such physics is undoubtedly present in the pnictide a^^'^^ , 
quantifying it would require a first-principles calculation 
of the relevant parameters. However, as our goal is to 
propose a simple mechanism that explains both the SPT 
and the resistivity anomaly, we focus on a more easily 
quantifiable approach to orbital ordering based instead 
on the Coulomb repulsion. Indeed, what our work in- 
dicates is that there is a rich set of models which can 
lead to orbital ordering in the pnictides. Our model 
is sufficiently simple and general that warrant its being 
taken seriously. The key insight gained from this study is 
not the detailed microscopic mechanism for this orbital- 
ordcring-induced SPT, which is rather standard^^, but 
its direct consequence - a resistivity anomaly, which can 
be captured by our model in quantitative agreement with 
the experimental results [see Fig.[3Ub)]. Furthermore, the 
stripe-like SDW and recently discovered anisotropy!^ of 
the magnetic exchanges naturally arise in our theory. 



II. ORBITAL ORDERING 

As being emphasized by pioneering earlier worbi^, the 
orbital degrees of freedom are important in the iron pnic- 
tides, which are intrinsically multi-orbital systems. For 
the Fe atom located at the center of the tetrahedron 
of four neighboring As atoms, the five d orbitals are 
spht into two groups, t2g (d^y, dxz, dyz) and Cg (d^2_j^2, 
d^^-2_^-2). Three of the five orbitals, dxy, dx2-y2 and 
d^z^-r^ are rotationally symmetric in the xy plane. So 
they are unlikely to have any effect on the SPT which is 
asymmetric in the xy plane. Then the only two possible 
candidates are the dxz and dyz orbitals. We propose the 
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FIG. 1: (Color online) (a) Equal numbers of dxz and dyz 
with a square lattice configuration, (b) Entirely dyz state 
with a < b. (c) Entirely dxz state with a > b. d) Relative 
energy difference A between configurations (a) and (b), or 
(c), as a function of lattice distortion 5. (Inset: an Ising- 
type transition where the order parameter M is defined as 
the difference between the numbers of occupied dxz and dyz 
orbitals.) 



following mechanism for the SPT, assuming these two 
orbitals are localized. At high temperature T > Tspt, 
dxz and dyz orbitals are degenerate, with equal numbers 
of electrons on both. A possible configuration is shown 
in Fig. [H^a), in which a square lattice is preferred. At 
low temperature, T < Tspt, there is a majority of ei- 
ther 

dxz 01" dyz- For dyz orbitals, the Coulomb repulsion 
of two neighboring sites is stronger along the y-direction 
than along the x-direction, which leads to a rectangular 
lattice with a < 6 as shown in Fig. [ijb), where a and 
b are unit lengths in the x and y direction, respectively. 
Similarly, when dxz dominates, the system will form the 
configuration of Fig. [IJc), which is degenerate with (b) 
by a rotation of 90 degrees. 

To demonstrate the viability of this mechanism, we 
need to compare the energies of configuration (a) and 
(b) in Fig. [TJ For simplicity, only the nearest-neighbor 
Coulomb repulsions arc considered. 
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\^^{r-R,)\^\i,p{r' ~R,f (1) 



where V'q('" — Ri) is the wave function of the a {a — 
dxz, dyz) orbital electron at site Ri. This integral can 
be evaluated by an importance-sampling Monte-Carlo 
method. In configuration (a), we choose a — b = qq ^ 
2.85 A, which is the typical experimental valuei for the 
1111-family. For configuration (b), we define the lattice 
distortion S as a = uq — S, and set b = Uq/q to keep the 
area of a unit cell constant. We calculate the relative 
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as a function of 5, where Ua and Ub are energies of config- 
urations (a) and (b) respectively. The results are shown 
in Fig. [T]^d). For a lattice distortion 0.03 k< 5 < 0.09 
A, the rectangular lattice (b) or (c) is more energeti- 
cally favorable. It is noted that this value is larger than 
the experimentally observed distortion of about 0.01 A^. 
However, the localized states are probably neither dxz 
nor dyz , but some combinations of the d orbitals, or even 
involve hybridization with As p orbitals^. Thus the pre- 
cise value of the distortion length can be smaller by tak- 
ing these factors into account. As already mentioned, 
other possibilities may also induce this ferro-orbital or- 
dering and the subsequent SPT. For example, Kriiger et. 
aft^"* derive a Kugel-Khomskii spin-orbital model and the 
resultant phase diagram does contain the same orbital 
configuration as proposed here. However, constructing 
the complete microscopic Hamiltonian that incorporates 
all the important physical processes requires a detailed 
knowledge of the relevant coupling parameters, which is 
currently unavailable. Thus the key point of our study 
is to put forth a simplified picture based on the coupling 
only to the Coulomb interaction in which the rectangular 
lattice with ferro-orbital ordering emerges spontaneously 
at low temperature because of its lower energy. 

Our model allows us to make the following conclusion. 
Upon the onset of the phase transition, a lattice distor- 
tion breaks the degeneracy between dxz and dyz- By 
occupying either one of these two orbitals, the system 
forms a ferro-orbital-ordered state and lowers its energy. 
It is this orbital-driven Jahn- Teller effect that induces the 
SPT. Defining Mi — ±1 for site i occupied by dxz and 
dyz orbitals respectively, we can write down an effective 
Ising-typc Hamiltonian for the SPT 



SPT 



(3) 



where Jspt should be on the order of the transition tem- 
perature, Tspt- So the SPT belongs to the Ising univer- 
sality class, as shown in the inset of Fig. W[d), where the 
order parameter M is defined as M = Mi/N. 

Recently, angle-resolved photoemission experiments 
using a linear-polarized laser bearoi^ show that at low 
temperature, the Fermi surface at the Brillouin zone 
center is dominated by a single dxz or dyz orbital, de- 
pending on the distortions. In their subsequent local- 
density approximation (LDA) calculationsiS., it is found 
that the density of states of the dyz orbitals with a lat- 
tice configuration of a < 5 displays a peak around 0.5 eV 
from the chemical potential, which is just the localized 
state predicted in our SPT model. A recent optical 
measurementiS, also suggests evidences of the orbital or- 
dering. 
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FIG. 2: (a) A schematic of a double- well potential, as the clas- 
sical analog of the two-level system, (b) Two types of scatter- 
ing processes between the itinerant electrons and the localized 
states. V^: diagonal scattering; and V'^: off-diagonal scatter- 
ing. 



III. RESISTIVITY ANOMALY 

The ferro-orbital-ordering-driven SPT mechanism has 
an important consequence, namely the resistivity 
anomaly. The essential physics is that of a Kondo prob- 
lem. The scattering of the itinerant electrons off two 
otherwise degenerate orbitals, dxz and dyz, will be sup- 
pressed by the gap opening, which results in a sharp drop 
of the resistivity upon the onset of the SPT. 

Above, TgpT, the two d^z and dyz orbitals are degen- 
erate. Below TsPT, the occupancy of the electrons in 
dxz and dyz orbitals becomes unbalanced as a result of 
the distortion of the crystal to configuration (c) [or (b)] in 
Fig. HI Thus, the electrons that remain in the dyz (or dxz) 
orbitals will have a higher energy and hence can lower 
their energy by jumping onto dxz (or dyz) orbitals. This 
process can be described by a localized two-level system. 
The classical analog, namely a double-well potential, is 
shown in Fig. [DJa). The corresponding Hamiltonian is 
given by 

+ J^o^alf^S/3a/3 (4) 



a/3 



where aj^ (cq) creates (annihilates) an electron on or- 
bital a and cr^^ is a Pauli matrix. We will choose an 
appropriate fictitious energy Aps to prevent the system 
from double occupancy. A is the energy splitting be- 
tween the two levels and Aq is the tunneling rate, as 
shown in Fig. H^a). By a rotation of the spin axis, this 
system can be diagonalized and the gap between the two 
eigenstates is E = ^/Aq 4- A^. 

As the parent compounds are actually metallic, there 
should be itinerant electrons present besides these lo- 
calized states. These two can be coupled as in 
the framework of the localized-itinerant dichotomous 



model o-^^'^°i^^'^^ . The starting Hamiltonian is^ 

H = H, + Htls + V (5) 
= ^i^fec^cfc. (6) 

i /ci<Tl,fc20-2 q/3 

where He, i?TLS and V represent the Hamiltonians for 
the itinerant electrons, the single two-level system and 
the interactions between the two, respectively. There are 
two kinds of scattering processes as shown in Fig. ^h) . 
One is the diagonal scattering described by the term, 
where the localized state remains on the same level. The 
other is the off-diagonal scattering initiated by the V 
term, where the localized state jumps onto the other 
level, is in fact zero, as it breaks time-reversal sym- 
metry. However, it should be noted that = does not 
hold for the renormalized vertex, since higher order terms 
are not necessarily local. We will also assume ^ 
as proposed previously^^. 

In fact, this system is very similar to the Kondo model, 
with the two orbitals dxz and dy^ representing the up and 
down-spin states on the magnetic impurity. We are go- 
ing to perform a similar scaling analysis following Ref . |23| . 
We define the dimcnsionless couphngs vl^f.^ — Vji^i^^No 
where A^o is the density of states at the Fermi level. Re- 
ducing the bandwidth from Dq to D and evaluating the 
vertex corrections up to the leading order, we have the 
scaling equations 



dv^fsju) 
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where are defined as v 



krk2 = J2flikl)Va0ff3{k2), 

with fa{k) being a complete set of spherical harmonics, 

faik) = i'-Y^{9k,(j>k), e'^" is the Levi-Civita symbol and 
u = ln(D/Do). We can express using the Pauli ma- 
trices as v^^p = v^al^p. Then the above scaling equations 
will be reduced to a set of coupled equations involving 
v^, yy and . These equations can be solved by sep- 
arating u into two regimes: (a) < <C and (b) 
yV ^ < w^. In regime (a), the solutions are 



v^(u) = w^(0)cosh[4w^(0)u] 
yy(u) = w^(0)sinh[4w^(0)u] 
v'-{u) ^ v'-{Q). 

In regime (b), we have 

[v^{u)f-[v%u)f^vl 

where vq is scale invariant and u^(u) satisfies 
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FIG. 3: (Color online) (a) Scaling of the coupling constants v' 
with respect to bandwidth D. (b) Energy gap and resistivity 
as a function of temperature T. (The experimental data of 
resistivity are extracted from Ref.^.) Setting the resistivity at 
T — 150K of our model equal to that of the experiment was 
the only fitting parameter. 



with the Kondo temperature identified as 



kBTk = Do 



w^(0) 



4w^(0) 



1/4-0^(0) 



(14) 



Using the parameters w^(0) = 0.33, v''{0)/v^O) = 
0.001, vy{0) = and Dq = 665 K^^, we obtained the scal- 
ing flows of w^, and shown in Fig.[3Ja), for £' = K. 
The corresponding Kondo temperature is Tk = 1.24 K. 
Reducing the bandwidth D, the system goes from weak 
to strong coupling. The resistivity due to the scattering 
of the two-level system can be calculated based on these 
renormalized vertices as in Ref.lU. At high temperature, 
we have two degenerate levels, d^z and dyz. When the 
temperature is reduced, the scattering from the states 
closer to the chemical potential increases, leading to a 
resistivity upturn of log T— as in the Kondo model. How- 
ever, upon the onset of the SPT, a gap opens between the 
two levels. If the bandwidth D is less than the gap E, the 
off-diagonal scattering is not allowed, since there are no 
states for the electrons to be scattered into. As a conse- 
quence, the scaling terminates at D — E. The electrons 
within the bandwidth E will no longer contribute to the 
resistivity. This is the mechanism behind the resistivity 
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FIG. 4: (Color online) Different possible spin configurations 
on a distorted lattice with a < b, which corresponds to the 
case that dyz is the majority orbital. 



anomaly. Our result is shown in Fig. [3l^b), which is in 
good qualitative agreement with experiment. We set the 
tunneling rate Aq = 2K, and the energy splitting takes 
the form 



Am^A(0)Jl-(^)' 



(15) 



where A(0) = Tspt = 150 K when T < Tspt- It should 
be noted that the overall behavior of the scaling flows and 
the resistivity are independent of the chosen parameters. 
This represents the explanation of the RA for the iron 
pnictides. 



IV. ORBITAL DRIVEN MAGNETISM 

Our model also offers a natural solution to the observed 
stripe-like anti-ferromagnetism. Before the SPT, we have 
an orbitally disordered state, in which the neighboring 
sites are occupied probabalistically by different orbitals. 
The resultant lack of overlap gives rise to a vanishing 
of any anti-ferromagnetic spin exchange and as a conse- 
quence no spin order. After the SPT, either dxz or dyz 
orbitals will dominate. Without loss of generality, we 
suppose that most sites are occupied by dyz , as shown in 
Fig. [TJb) . Due to the larger overlap of the wave functions 
on neighboring sites in the y-direction than that in the 
x-direction, the hopping integral ti, should be larger than 
ta- For the nearest-neighbor spin exchange, Ji ~ f^/U, 
we have that Jia < Jib- So the spins on the longer axis 
have a stronger tendency to be aligned oppositely. The 
spin configuration AFM2(b) in Fig. |4]is not favored. As 
has been suggested^'^^, we can further introduce a next- 
nearest-neighbor exchange J2- If J2 > •/ia/2, which is 
very likely for a relatively small JiJ^, AFM2(a) will 
have a lower energy than AFMl, as shown in Fig.lH and 
emerge as the ground state at low temperature, which 
has already been confirmed by the experiments^. In con- 
trast with other theories in which the SPT is induced by 
the spin degrees of freedoms, on this account, the forma- 
tion of the SDW is actually a result of the ferro-orbital 
ordering accompanying the SPT. 
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In fact, we are able to construct a universal Hamilto- 
nian describing both the SPT and SDW, with a spin-orbit 
coupling modeli^i'i^'i^, 

Hso = JsptJ2m,M,+ J2{M,,Mj)S,-Sj 

i 

+y) ' ^i-\-y 

(16) 

i 

where the spin exchanges are given by 

Ji^{Mi,Mj) = 5M,.M,{Juhu.i + JiJivh,-i) (17) 

Jly{Mi,Mj) = 5m,,M, {Jla^Mi,! + Jlb5Mi,-l) {^^) 

J2{M,,M,) = 5m,,m,J2 (19) 

where Mi, representing the orbital degrees of freedom, is 
defined to be ±1 for dxz and dyz respectively, as in SecHIl 
Clearly, in this model, the spin order will not occur until 
the formation of the ferro-orbital ordering at Tspt, which 
is on the order of Jspt- Below Tspt, the spin degrees 
of freedom can be described by an anisotropic Heisen- 
berg model, whose transition temperature to the spin- 
ordered state, Ts, would depend on the spin exchanges, 
Jia, Jib and J2. If Ts < Tspt, we would have two sepa- 
rate second-order transitions, Tsdw = Ts < Tspt, as in 
the case of the 1111-family. For the 122-family, which has 
a shorter Fe-Fe bond length, it is expected this would en- 
hance the spin exchange J, likely leading to Ts > Tspt- 
But the SDW will not form before the SPT, since there 
is no spin exchange until the SPT obtains. So there is 
only one first-order transition, Tsdw = Tspt- 

Furthermore, this anisotropic Heisenberg model has 
also been proposed on experimental groundsi^ to fit the 
spin-wave spectrum seen in the inelastic neutron scatter- 
ing data. Our theory gives a direct explanation for the 
observed anisotropy of magnetic exchanges. Note their 
resultsi^ do rely on a negative J", which is not obtained 
by our simple model. However, this difficulty can over- 
come by introducing a Hund's coupling between these 
localized spins and itinerant electrons. 



where cr^j^/ are the Pauli matrices. The hopping of the 
itinerant electrons with this Hund's coupling will give 
rise to an effective ferromagnetic coupling22i22ii2S between 
neighboring spins. After taking this into account, we 
will eventually have the spins on the shorter axis coupled 
ferromagnetically. The full details of this model arc the 
subject of a future study. 



V. FINAL REMARKS 

To conclude, we have proposed that the SPT and RA 
in the iron pnictides are due to the opening of a gap 
between two otherwise degenerate orbitals. While our 
mechanism for the structural phase transition is a stan- 
dard Jahn- Teller distortion driven by a minimization of 
the Coulomb repulsion, the key point of this paper is 
that the resulting simple two-level system can resolve the 
previously unexplained resistivity anomaly. The mecha- 
nism proposed here is independent of an applied mag- 
netic field as is seen experimentally^. Only in a fcrro- 
orbital-ordered state after the SPT does the stripe-like 
SDW form. This is the reason why these three phenom- 
ena, SPT, RA and SDW, are closely related and almost 
always coincide with one another. In doped materials, 
extra electrons or holes will break the uneven occupa- 
tions of dxz and dyz, thus diminishing the Jahn- Teller 
effect. So the SPT, RA and SDW wiU aU become less 
pronounced and shift to lower temperature, eventually 
vanishing at some critical doping. These are all observed 
experimentally, lending credence to our model. 

After this work was posted on arXiv, several similar 
paper o^"'^^ appeared, based on the same orbital physics 
we utilized here, which supports our theory that the or- 
bital ordering is the driving mechanism for the SPT, RA 
and SDW. 
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